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Abstract. The Circinus Galaxy has been observed, for 
the first time above 10 keV, by BeppoSAX. An excess 
emission above the extrapolation of the best fit 0.1-10 keV 
spectrum is apparent in the PDS data. The most likely 
explanation is that we are observing the nucleus through 
a column density of ~ 4 x 10 24 cm -2 , i.e. large enough 
to completely block the radiation below 10 keV but small 
enough to permit partial transmission above this energy. 

Key words: X-rays: galaxies - Galaxies: Seyfert - Galax- 
ies: individual: Circinus 



1. Introduction 

The Circinus Galaxy hosts one of the nearest Active Nu- 
clei (d ^4 Mpc) and has been therefore much studied after 
its discovery (Freeman et al. 1977), despite the low Galac- 
tic latitude and therefore the relatively large interstellar 
absorption (N H ,Gai ~3x 10 21 cm -2 ). 

The intense coronal lines (Oliva et al. 1994, Moor- 
wood et al. 1996), the [O III] ionization cone (Marconi 
et al. 1994), the variable H2 maser emission (Greenhill et 
al. 1997), the broad Ha line in polarized light (Oliva et 
al. 1998a), and the reflection-dominated 2-10 keV X-ray 
spectrum (Matt et al. 1996, hereinafter M96) support the 
Seyfert 2 classification for the Circinus Galaxy. The in- 
tegrated IR luminosity is about 4xl0 43 erg s -1 , a large 
fraction of it being probably reprocessing of the nuclear 
radiation (Moorwood et al. 1996; Maiolino et al. 1998). 



The ASCA observation (M96) revealed a flat X-ray 
spectrum with several emission lines, including a very 
prominent iron Ka line at about 6.4 keV. These findings 
were interpreted by M96 in terms of reflection by cold 
matter (possibly the inner wall of the putative torus en- 
visaged by unification models, e.g. Antonucci 1993) of an 
otherwise invisible nuclear emission. The absence of any 
transmitted component in the ASCA band (i.e. up to 10 
keV) indicates a column density for the absorbing matter 
in excess of 10 24 cm -2 , i.e. a Compton-thick source. 

The Circinus Galaxy was observed by BeppoSAX 
(Boella et al. 1997) on 1998 March 24 as part of a 
Core Program aiming to study the broad band spectra of 
Compton-thick Seyfert 2 galaxies. Here we report on the 
spectral analysis of the source, with particular emphasis 
on the high energy part of the spectrum. 

2. Observation and data reduction 

We will discuss here data from three instruments: the 
LECS, the MECS and the PDS. The MECS is presently 
composed of two units (after the failure of a third one), 
working in the 1-10 keV energy range. At 6 keV, the 
energy resolution is ~8% and the angular resolution is 
^0.7 arcmin (FWHM). The LECS has characteristics sim- 
ilar to that of the MECS in the overlapping band, but its 
energy band extends down to 0.1 keV. The PDS is a pas- 
sively collimated detector (about 1.5x1.5 degrees f.o.v.) 
working in the 13-200 keV energy range. The effective ex- 
posure time of the observation was 1.37xl0 5 sec (MECS), 
8.37xl0 4 sec (LECS) and 6.32 x!0 4 sec (PDS). 
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A second source is clearly visible in the LECS and 
MECS images at about 5' from the Circinus Galaxy. This 
source was discovered by ASCA (M96), and is also present 
in the ROSAT HRI image (Ward 1998). The nominal po- 
sitions of both sources, as measured by BeppoSAX, are 
shifted by the same amount (about 2') and in the same 
direction with respect to the HRI positions (which in turn 
are consistent with the optical position, for Circinus, and 
the ASCA position, for the second source). Only the Y 
and X startrackers, one gyro and the solar sensor were in 
use during the observation, but not the Z startracker, i.e. 
that coaligned with the telescopes. In this configuration, 
absolute position reconstruction can be in error up to 3'. 
The second source is about 3.7 times fainter in both in- 
struments. The spectrum may be well fitted either by a 
power law with T ^2.35 and Nh ~ 1-3 x 10 22 cm -2 or 
a plasma model with kT ~6 keV and Nh ~ 3.4 x 10 21 
cm -2 . In both cases, the extrapolation to the PDS band 
lies more than two orders of magnitude below the data, 
and its contribution to the PDS data may be neglected. 
Of course, in principle the excess emission in the PDS 
band, described in Sec. 3, could be due to this source, if it 
is heavily absorbed (Nh — 10 24 cm -2 ). This would imply 
a physical scenario similar, but much more extreme and 
then less plausible, to that discussed below for Circinus. 

To avoid, as far as possible, contamination from the 
second source we extracted data for both the LECS and 
the MECS in 2' radii centred on the source. The back- 
ground subtraction was performed using blank sky spec- 
tra from the same region of the detector field of view. 
The resulting count rates are 4.41 (±0.08) x 10 -2 in the 
LECS and 1.050(±0.009) x 10" 1 counts/s in the MECS. 
PDS data have been reduced following the standard proce- 
dures described in Matt et al. (1997). In addition, we em- 
ployed crystal temperature dependent Rise Time thresh- 
olds, which lead to a reduction by a factor up to 50% 
of the instrumental background. The PDS count rate is 
2.01±0.04 counts/s. 

The spectra were fitted using the calibration matri- 
ces released on September 1997. In the fit procedure, the 
LECS normalization was left free to vary with respect to 
the MECS one, while the PDS normalization was fixed to 
0.8 (Grandi et al. 1997). The residual uncertainty on the 
PDS absolute flux calibration (~ ±10%) does not affect 
substantially the following results. All quoted fluxes refer 
to the MECS. Errors correspond to 90% confidence level 
for two interesting parameters (i.e. Ax 2 =4.61). 

Meaningful LECS and MECS light curves of the Circi- 
nus Galaxy cannot be easily constructed: due to the prob- 
lems in the aspect reconstruction, a small extraction ra- 
dius results in a fictitious variation of the count rate due to 
the wobbling; if a large extraction radius is instead chosen, 
the contamination from the second source is significant. 
In the latter case, a variation up to 10% on time scales of 
hours is observed in the MECS. If this is fully attributed 
to the second source, as should be if the interpretation 



of the Circinus spectrum below 10 keV as pure reflection 
component is correct, this would imply a ^50% amplitude 
variations in this source. Finally, the PDS light curve is 
consistent with being constant (with a ^40% upper limit 
to the variability amplitude). 

3. Spectral analysis and results 

We firstly analysed the spectrum below 10 keV, to verify 
the consistency with the ASCA results (M96). The spec- 
trum is very flat with a huge iron Ka line, as observed 
by ASCA and interpreted by M96 as reflection dominated 
emission, the reflection probably occurring in circumnu- 
clear matter (hereinafter identified, for simplicity, with the 
"torus" familiar from the unification models). The 2-10 
keV flux is 1.4xl0~ n erg cm~ 2 s -1 , consistent with the 
ASCA one. At least other five lines, to be identified with 
Si, S, Ar and Ni Ka and Fe K/3 lines, are required, as well 
as excess emission below a few keV. The properties of the 
emission lines, as well as of the soft excess, are discussed 
in detail elsewhere (Guainazzi et al., 1998). 

When the model described above is extrapolated to the 
PDS band, a huge excess emission is apparent (see Fig.0). 
No known bright sources are present in the field of view 
of the PDS. The Galactic longitude (311°) and latitude 
(-3.8°) of the source are about 10° outside the bulk of the 
Galactic Ridge emission (Warwick et al. 1985). Following 
the modelling of Yamasaki et al. (1997), we found that 
the contribution of the Galactic Ridge to the PDS count 
rate is at most 7%. Therefore, it is very likely that we 
are observing a further component in the Circinus spec- 
trum, possibly the nuclear radiation piercing throughout 
a very large column density absorber, as in the case of 
NGC 4945 (Iwasawa et al. 1993; Done et al. 1996) and 
Mkn 3 (Cappi et al. 1998) In order to conceal itself below 
10 keV, the absorber column density must exceed 10 24 
cm -2 , i.e. must be Compton-thick. This implies that the 
photons which eventually emerge may have suffered one 
or more Compton scatterings. We have therefore adopted 
a model, based on MonteCarlo simulations, which prop- 
erly treats photon transfer in very high column density 
matter (work in preparation), including Compton down- 
scattering and the Klein-Nishina decline. We assumed a 
spherical geometry as the simplest approximation of the 
(actually unknown) distribution of the matter. If a plane- 
parallel geometry, with normal illumination, is employed, 
the results do not change significantly. 

We have then fitted the whole spectrum with the fol- 
lowing model: 

F(E) = [A 1 Tr(E,N Utl ,r h ,E c )+A 2 R(E,r h ,E c ) 

+AzE- r ° + emission lmes]e~ Nli - 2a ^ (1) 

where 

Tr= [ G(E,E ,N u>1 )Eo r *e(- E °/ E JdE (2) 

J E 
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1. LECS+MECS best fit model extrapolated to the PDS 
The model is composed by a pure reflection component 
1.60), a soft power law (F a = 1.45) plue the emission lines. 



Moorwood et al. (1996) and Oliva et al. (1998b) based on 
the IR lines. The X-ray to OIII ratio is typical for Seyfert 
galaxies (Bassani et al. 1998). 
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and G{E,Eq) is the Green function relating input and 
transmitted spectra. Th and T s are the power law indices 
of the hard (nuclear) component and of the soft compo- 
nent, respectively; -/Vh.i is the column density of the torus 
and JVh,2 that of any external matter along the line of 
sight. R(E,Th, E c ) is the pure reflection spectrum, ob- 
tained using PEXRAV in XSPEC fixing cos(#) to 0.45 and 
rel_ref 1 to -1. 

The 20-100 (50-150) keV total flux is 2.3 xl0~ 10 
(0.9xl0~ 10 ) erg cm -2 s _1 , which makes the Circinus 
Galaxy one of the brightest AGN in hard X-rays. The 
best fit parameters are summarized in Table 1 (model 1), 
while the best fit model is shown in Fig. |[ T s and Th are 
very similar, and consistent within the errors. It is possible 
that the soft emission is due to scattering of the nuclear 
radiation by ionized matter (Guainazzi et al. 1998). 

The normalization of the nuclear component, A±, is 
about 5 times larger than that of the reflection component, 
Ai\ therefore, the solid angle subtended by the reflecting 
matter to the nucleus and visible at infinity is ~0.2x2-7r. 
Assuming a torus geometry like that adopted by Ghisellini 
et al. (1994) (i.e. a conical geometry with half-opening 
angle of 30°), we derive an inclination angle ^ 40°. 

The 2-10 keV nuclear luminosity, after correction for 
absorption, ranges from 3.4 xlO 41 to 1.7xl0 42 erg s _1 (as- 
suming a distance of 4 Mpc). The IR luminosity (as de- 
fined by Mulchaey et al. 1994) is 20 to 100 times larger. 
This IR/ (2-10 keV) ratio is somewhat higher than the av- 
erage value for Seyfert galaxies, but by no means unique 
(see Fig. 2a of Mulchaey et al. 1994; Awaki 1997). If a sig- 
nificant fraction of the IR emission is reprocessing of the 
nuclear component, a large UV nuclear emission is im- 
plied, qualitatively in agreement with the calculations of 
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Fig. 2. Best fit reflection + transmission model (in vF v ). 

We have also tested a different scenario (model 2 in 
Table 1), in which the transmitted component is simply a 
fraction of the total reflection component, i.e. the reflect- 
ing matter is partially covered. This case would represent 
a scenario in which the line of sight towards the nucleus is 
completely blocked (i.e A^h in excess of ~10 25 cm -2 ), but 
the optical depth diminishes at higher torus latitudes, so 
permitting transmission of some of the radiation reflected 
by the far side of the torus. The spectrum is similar to 
that of eq. 1, but substituting, in the integral of eq. 2, the 
reflection component to the power law one. From the sta- 
tistical point of view, this model is as good as model 1. The 
A\ I A2 ratio is now the ratio between the solid angles sub- 
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tended by the covered and uncovered matter to the source 
and visible at infinity. From Table 1, it can be seen that 
the covered fraction of the reflector is 5 times larger than 
the uncovered one. No estimate of the nuclear luminosity 
is possible in this model, and therefore no information on 
the torus geometry can be derived. 

The observed value of the torus column density is 
4.3xl0 24 cm -2 in model 1 (and not much lower in model 
2), i.e. about 350 times larger than the value expected 
from the extinction in the V band (Ay =5 mags) measured 
from narrow lines Balmer decrement (Oliva et al. 1994) 
assuming a Galactic dust-to-gas ratio. Very likely, the 
absorbers responsible for the obscuration of the nucleus 
and the narrow line region are different, the former being 
distributed on much smaller scales. Granato et al. (1997) 
fitted the ISO spectrum (Moorwood et al. 1996) with their 
torus model, and found Ay =50 mags for the torus, corre- 
sponding to a Ah more than one order of magnitude lower 
than our measured value, thus implying a dust depleted 
material. If their estimate is correct, the absorbing matter 
should then be located inside the dust sublimation radius 
which, given the nuclear luminosity derived from model 
1, is of order of 0.01 pc, very close indeed to the central 
engine and suggestive of a possible identification of the 
absorbing material with the BLR. 

Finally, a high energy cutoff with e-folding energy be- 
low 100 keV is required by the data in both models. This, 
together with the rather flat intrinsic power law, make 
Circinus similar, as far as the primary spectral shape is 
concerned, to NGC 4151 (see e.g. Piro et al. 1998). 

4. Discussion 

The main results of the BeppoSAX observation of the 
Circinus Galaxy can be summarized as follows: 

i) The spectrum below 10 keV is dominated by a reflection 
component, in agreement with the ASCA result (M96). 

n) The extrapolation of the LECS+MECS best fit spec- 
trum falls short of the PDS data. A further component 
transmitted through a ^4xl0 24 cm~ 2 absorber is re- 
quired, but it is not possible to establish, on statistical 
ground, whether we are looking at nuclear radiation or at a 
reflection component passing through the outer edge of the 
torus, even if the rather ad hoc geometry required in the 
second scenario makes the first hypothesis more appealing. 
To further address this issue, one has to search for short 
term variability, a task which requires much more sensitive 
hard X-ray detectors, like those on board Constellation- 
X. If the transmitted radiation is actually the nuclear one, 
by comparing the nuclear luminosity and that of the re- 
flection component the inclination angle of the system is 
constrained to be less than 40°, if the torus geometry pro- 
posed by Ghisellini et al. (1994) is adopted. 

Finally, it is worth noting that two of the nearest 
AGNs, the Circinus Galaxy and NGC 4945, both have 



absorbing matter with Ah between 10 24 and 10 25 cm~ 2 , 
suggesting that this kind of source is rather common in 
the Universe. This is particularly interesting because with 
such an absorber most of the energy in the X-ray band 
comes out at a few tens of keV (Fig. 2), where the energy 
density of the Cosmic X-ray Background peaks. 
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